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Abstract 



We show that at certain values of QCD condensates the nucleon QCD sum rules 
with "pole+continuum" model for the hadron spectrum obtain an unphysical solution. 
This provides constrains for the values of condensates to be consistent with existence 
of a physical solutions. The constrains become much weaker if the radiative corrections 
are included perturbatively. We demonstrate that the most important dependence of 
nucleon mass on the quark scalar condensate becomes less pronounced under factorization 
assumption for the four-quark and six-quark condensates. 

1 Introduction 

The QCD sum rules invented by Shifman et al. [1] enable to express vacuum characteristics 
of hadrons in terms of the vacuum expectation values of QCD operators. This approach was 
employed to description of nucleons in [2] and [3, 4]. The improved analysis was presented later 
in [5, 6]. The method was used also for description of delta- isobars [5, 6] and of the baryons 
containing heavier quarks [7, 8]. Further applications of the vacuum nucleon QCD sum rules 
are reviewed in [9]. The method was expanded also for the cases of finite temperatures [10] and 
densities [11]. 

The main tool of the vacuum QCD sum rules for a hadron is the dispersion relation for 
polarization operator 



with j{x) a local operator ("current") carrying the quantum numbers of the hadron. It is 
nucleon (proton) in our case. The dispersion relation is considered at large values of \q^\{q^ < 0) 
where n(g^) can be represented as a power series in with the vacuum expectation values of 
QCD operators as the coefficients of the expansion. This is known as operator power expansion 



Recall the main milestones of the QCD sum rules analysis. Following [2] one can write 
dispersion relations 




(1) 



(OPE) [12]. 




(2) 



{i = q, 1) for the ingredients n*(g) of the polarization operator 

U{q) = g^7''n<^ + /n^(g2). (3) 

We equal the OPE of the left-hand side (LHS) of Eq. (2) to the contribution of the observable 
hadrons to its right-hand side (RHS). The latter is usually approximated by the "pole + contin- 
uum" model in which the lowest pole is written exactly while the other states are approximated 
by continuum: 

- Im n^(P) = Xl6{k^ - m^) + e{e - W^)r{e) . (4) 
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Of course, the "pole + continuum" model is reasonable only if the contribution of the pole 
exceeds that of the continuum. 

Behavior of both sides of Eq. (2) prompts the choice [1] 

f{e) = - Im (OPE)(^2) 

Note that in this approach the continuum threshold W'^ does not coincide with the physical 
continuum threshold. Thus the position of the lowest pole m, its residue and the model 
continuum threshold W'^ are the unknowns which are expected to be determined by the QCD 
sum rules equations. The standard next step is the Borel transform, after which Eqs.(2) take 
the form 

£«(M2) = i?«(M2); C'iM^) = R\M^). (5) 

Here i2*(i?*) are the Borel transforms of the LHS (RHS) of Eq. (2), is the Borel mass. This 
approach provided good results for the nucleon mass and for the other nucleon parameters [9] . 

Note that both D{M'^) and K^^M^) are calculated in framework of certain models. The 
OPE expansion for 0{M'^) is increasingly true at large values of M^. The "pole -|- continuum" 
for model for i?*(M^)is increasingly true at small values of M^. Important assumption is that 
there is a region of intermediate values of where both approximations work and reproduce 
to some extend the true (unknown) function of M^. Thus our task is to find the interval of 
the values of M^, where the functions £*(M^) can be approximated by the functions i?*(M^) 
and to find the set of parameters m, and W"^ which insure the most accurate approximation 
of D{M'^) by i?*(M^). The set of values of parameters m,\\ and W"^, which minimize the 
function 




X\m, A^, ly^) = E E ^ V,v..2^ ' (6) 



will be referred to as a sohition of the sum rules equations. 

Note that it is important to obtain "duality" between the LHS and RHS of Eq.(5) in some 
interval of the values of M^, but not at certain point M?. Therefore we will look for the three 
unknown parameters simultaneously. 

Both the interval of the values of the Borel mass (" Borel window" ) and the solution of the 
sum rules depend on the form of the proton nucleon current j (a;), which is not determined in 
an unique way. The general form is [2, 13] 

j{x]t) = ji{x) +tj2{x), (7) 
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with 

jl{x) = eabc[Uaix)Cdb{x)]-f5Uc{x); j2{x) = eabc[Uaix)Cj5db{x)]Uc{x), 

where u and d are the quark operators, a, b, c are the color indices, T denotes a transpose and 
C is the charge conjugation matrix, while t is an arbitrary parameter. 

Following [2], we shall use the current determined by Eq.(7) with t — —1. It can be written 
(up to a factor 1/2 ) as [2] 

= eabc[Uaix)Clf,Ubix)]'j5Ydcix) (8) 

This choice was shown in [14] to be most relevant for description of nucleons, since the polariza- 
tion operator Il{q^) calculated with this current satisfies two main requirements. On the RHS 
of Eq. (5) the contribution of the nucleon pole exceeds that of the higher states (approximated 
by continuum). On the LHS of Eq.(5) the higher order terms of series should drop fast 
enough to be consistent with the convergence of the OPE. It is important also that there is a 
gap between the position of the lowest pole and the effective threshold W"^. This choice of 
the current was advocated in [13]. There are many papers, in which this very current was used 
- see, e.g., [15]. 

Several lowest terms of the OPE for the current (8) have been calculated in [2, 5]. The 
leading term depends on q"^ as g^lng^. It comes from the free three-quark loop. The higher 
order OPE terms contain the matrix elements 

{0\qq\0), {0\-Gl,GlM, mqqq\0), 

TT 

etc. Analysis carried out in [5] contained also the most important radiative corrections in which 
the QCD coupling constant is multiplied by "large logarithm" Ing^. Corrections of the order 
(a^lng^)" to the leading OPE terms have been calculated earlier in [16]. 

The appropriate interval 

0.8 GeV^ < < 1.4 GeV^ (9) 

(" Borel window" ) was found in [5] . The values of the lowest OPE terms at conventional values 
of the QCD condensates enable to expect the convergence of OPE series. Also, for the solution 
found in [5] the contribution of the pole exceeds that of the continuum. 

The OPE series for 11^ and 11^ start from the terms q^ In and q"^ In q"^ correspondingly. 
Thus, in somewhat straightforward interpretation of OPE only 11'^ contains a leading term. 
This allows to consider the chirality conserving structure 11^ as a more important one. However 
actually we consider the values of the Borel mass of the order of the proton mass. After the 
Borel transform the leading contributions to 11^ and U.^ /m are of the same order. The chirality 
violating sum rule for 11^ can be considered as important as that for 11''. The two sum rules 
were considered on the same terms, requiring the same accuracy for both of them. That's why 
the terms corresponding to 11'^ and 11^ on the right hand side of Eq.(6) were included with the 
same weights. 

A weak point of this procedure is that the choice of parameters m, and W"^ may be not 
simple. There can be several local minima of corresponding to several sets of the parameters. 
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Note that the values of QCD condensates are known with large uncertainties. The expec- 
tation value (0|gg|0) can be determined with the larger accuracy than the other condensates 
due to the Gell-Mann-Oakes-Renner relation [17]. There is no experimental data on values of 
four- and six-quark condensates. They can be calculated in the factorization approximation 
for the case of large number of colors A'^c ^ 1. The accuracy of this approximation for = 3 
is obscure. The value of gluon condensate was initially obtained from the QCD sum rules 
for p mesons [18]. However the sum rules for other mesons lead to somewhat smaller [19] or 
larger [20] values, with the latest analysis presented in [21]. Hence it is reasonable to study 
dependence of the nucleon parameters on the values of the QCD condensates. 

This dependence should be investigated together with inclusion of radiative corrections. The 
latter imitate modification of the values of the condensates providing the contributions ~ 
and as In to each of OPE term. Thus, including the radiative corrections into our analysis, 
we can separate the effects of uncertainties in the QCD condensates values. 

In the present paper we analyze the role of the unphysical solution for the nucleon QCD 
sum rules, which corresponds to the continuum contribution exceeding that of the pole. We 
mentioned this solution in our earlier papers [22, 23]. Here we show that both physical and 
unphysical solutions provide local minima of the function ^"n^ W'^) determined by Eq. (6). 

For the absolute values of condensates, which differ noticeably from the conventional values 
(still consistent with convergence of the OPE series on the LHS of the sum rules), the minima 
corresponding to physical solutions may vanish. 

We show that inclusion of the radiative corrections modifies the situation. One could expect 
this, since it was shown in [22] that the radiative corrections affect mostly the value of the 
nucleon residue A^. After the corrections of the order as are included perturbatively, the 
physical solution exists for a broader interval of the values of condensates. Also, domination 
of continuum contribution over that of the pole for the unphysical solution becomes stronger. 
The unphysical solution becomes "more unphysical". 

Note that the problem of the pole dominance emerged in other QCD sum rules studies. In 
review on the QCD sum rules analysis of the pentaquark states [24] the authors found that 
it is very difficult to satisfy the requirements of pole domination, OPE series convergence and 
stability within the Borel window simultaneously. This contrasts the earlier statements (cited in 
[24]) that the QCD sum rules support the existence of the pentaquark. However, the authors of 
[24] do not make a definite statement on the QCD sum rules predictions about the pentaquark 
states. In the QCD sum rules analysis of the light tetraquark states [25] it was found that the 
solution with the domination of the pole can be obtained only for small values of the Borel 
mass M^, where the OPE series does not converge. The authors of [25] conclude the QCD sum 
rules analysis does not support existence of light tetraquark particles. In view of the analysis 
carried out in the present paper, inclusion of the radiative corrections may become important 
here. 

In Section 2 we analyze the interplay of the physical and unphysical solutions taking into 
account only the leading radiative corrections. We include corrections of the order ag in Sec. 3. 
We summarize in Sec. 4. 
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2 Interplay of the physical and unphysical solutions 



We start by representing the nucleon QCD sum rules [5] without using the factorization hy- 
pothesis for the condensates of the high dimensions. The LHS of the Borel transformed nucleon 
sum rules (Eq.(5)) with inclusion of the anomalous dimensions (i.e. of the corrections of the 
order {aslnq'^)"') can be written as 



(10) 



Here the lower indices show the dimensions of the condensates. If the current j in Eq.(l) is 
given by Eq.(8), the terms on the right hand sides of Eq. (10) are [2, 5] 



M^Eo 



4L 

C7 



^6 = g Ce-L , 



A ^---^ 
^ 3 M2 



12 ' " 81 M2 
Here L accounts for the leading radiative corrections ~ In [16] 



(11) 



dngVA^ 



QCD 



with the anomalous dimension ^ — 4/9 {L — lif the leading radiative corrections are neglected). 
After the Borel transform 



The condensates q are 

C3 = -(27r)^(0|^g|0), 

C6 = (27r)4(0|gggg|0), 



C4 = (27r)V0 



a. 



TT 



C7 = -(2ir)V0 



2, 



q—G''q\0 

TT 



Cs^{2nn0qqq-Gl,-a^.q 



TT 



C9 = -{2nf^ {0\qqqqqq\0) ■ 

TT 



(13) 



Note that the structure of the condensates cg and Cg is indeed more complicated. For ex- 
ample, Ce contains the condensates {0\qT^qqT^q\0) with being the basic 4x4 matrices 
(F^^ = I,^n,^5:i^5^ij,,crfj_^). The same is true for Cg. Hence the matrix elements {0\qqqq\0) 
and {0\qqqqqq\0) in expressions for cg and Cg in Eq. (13) are somewhat "effective" 4g and 6q 
condensates. Denote 



C3 

C7 = aobofqg , 



Cg — — (J'oJeq 

TT 



(14) 



Here we introduced dimensionless parameters /j. In the factorization approximation f^q = 
fqg = f2qfb, fsq = fiq- ^c shall investigate the dependence of nucleon parameters on QCD 
condensates modifying the values of /j. Note also that Cg = ulal with nl a; 0.8 GeV^ [2]. 
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Following [2] we write 



a = -(27rf(0|gg|0) , 6 = (27r)^(0 ^ 0) . (15) 

For traditional choice of the normalization point // = 0.5 GeV 

a = oo = 0.55 GeV^ 6 = 6o = 0.5GeV^ (16) 
Putting all /j = 1 in Eq. (14) we come to the standard nucleon sum rules [2, 5] with 

^0--^, ^4-^^, A,--aL, A,---— a, 

3 ^ ' ^ 12 ' ^ 81 TT M2 ^ ^ 



Here 



£;„ = £;„(x) = i-e-^5^ — , 71 = 0,1,2, x^wyM\ 

i.e. the contributions of continuum are transferred to the LHS of Eq. (5). The RHS of Eq. (5) 

i?^(M2) = e— ^ R\M^) = mA^e-"^'/^' (18) 

contain only the contribution of the nucleon pole with A^ = 32 7r'^A^. 
In the Borel window 

with = 2aM^ is just determined by Eq.(17) before the threshold contribution is trans- 
ferred to the LHS. (Note that A4 and B^ are numerically small due to a small coefficient 
connected with contributions of the gluon condensate). This is consistent with the hypothesis 
about convergence of the OPE series. 

We can expect that the convergence will not be spoiled by the disregarder terms. Note 
that the term can be viewed as coming from expansion of the expectation value of the 
operator q{0)q{x) in powers of x"^. The ratio As/Aq is the characteristic scale for further 
expansion in powers of M"^, caused by expansion in powers of x^. The mixed quark-gluon 
condensates of higher dimensions are expected to have small numerical factors, connected with 
the gluons, similar to C4 and Cy. Finally, there is the only QCD parameter Aqcd, a-nd there 
can be contributions containing the factor Ag^-^^/M^ ^ 1. 

In framework of the factorization hypothesis the method developed in [2, 5] provides 

m = 0.931 GeV , A^ = 1.86 GeV^ = 2.09 GeV^ , (20) 

if the numerical values (16) are employed. We assume Aqcd ~ 150 MeV. In the one-loop 
approximation this corresponds to 0^5(1 GeV^) ~ 0.37, which is consistent with the PDG data 
[26]. 
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If the radiative corrections are neglected, i.e. L = 1 the solution appears to be 

m = 0.930 GeV, = 1.79 GeV^ = 2.00 GeV^ (21) 

One can see that for the solution represented by Eq. (21) the contribution of the pole exceeds 
that of continuum more than twice. 

However the successive inclusion of the OPE terms leads to another solution. Taking into 
account only the condensates with the dimensions d — 3, 4 we find a trivial solution m — 

= = 0. Inclusion of the condensate with d = 6 keeps m = 0, — but provides 
A2 = 4/3 = 0.4 GeV^ The condensates with ci = 7 and ci = 8 require the nonzero values of 
m and W^, i.e. 

m = 0.6 GeV , A^ = 0.79 GeV^ , = 1.0 GeV^ . (22) 

We treat this solution as an unphysical one since in the Borel window determined by Eq. (9) the 
contribution of the pole is smaller than that of the continuum. For example, at the characteristic 
value = 1 GeV^ of the Borel window the ratio of the continuum and pole contributions is 
2 and 1.75 for the 11^ and 11^ structures correspondingly. 

As we shall see below, the unphysical solution manifests itself when the QCD condensates 
deviate from their conventional values. Note that for all the cases discussed below Eq.(19) is 
still true, and thus the convergence of the OPE series is not violated. Also (see, e.g. Fig.l 
below) the value of is small enough. Hence, the only reason for assuming this solution to be 
unphysical is the domination of the contribution of the continuum over that of the pole. 

It is instructive to trace the dependence of nucleon parameters on the value of gluon con- 
densate. At fb — 1, i.e. at the value of the gluon condensate corresponding to Eq. (16) the 
functional (6) has two local minima corresponding to the solutions (21) and (22). The deeper 
minimum is provided by the unphysical solution (22). At < 1 we still have two minima, and 
the values of the nucleon parameters do not change much even for fb = due to a relatively 
small value of the term A4 in Eq. (17). However for > 0.2 the deeper minimum corresponds to 
the unphysical solution (22). Somewhat straightforward employing of the chi-squired method 
may leave the physical solution unnoticed. The situation is more dramatic for /{, > 1. At 
fb > 1.04 the certain minimum corresponding to the physical solution vanishes and only the 
unphysical solution survives. 

In order to illustrate the role of the unphysical solution in the Borel window we introduce 
the function 

with jC^ and defined by Eq.(lO). Using Eq.(18) we see that for the values of determined 
by Eq.(9) we can expect m(M^) ^ const — m for the solutions of the sum rules equations. 
Putting fb = 0.6 (to make the difference of the corresponding values more visible) we expect 
the solutions to be close to those described by Eqs.(21),(22) for fb = 1. The functions m{M'^) 
for = 2.09 GeV2- see Eq.(21) and for = 1.0 GeV^ -see Eq.(22) are shown in Fig.3. One 
can see that the unphysical solution exhibits a more stable behavior then a physical one. 

In Fig. 3 we show the dependence of m, A^ and W"^ on the value of /;,. 
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Note that if there is an unphysical solution which corresponds to the absolute minimum of 
the functional x^(m, A^,l^^) defined by Eq.(6), the search for another solution and its inter- 
pretation becomes more complicated. It is rather simple if both minima are sharp. However, if 
the unphysical solution corresponds to a wide minimum, while the second solution corresponds 
to a shallow one, interpretation of the latter solution requires additional analysis. 

Now we come to variation of the values of the quark condensates. Let us first modify 
the value of the condensate (0|gg|0), keeping the other ones to be unchanged. Hence we put 
fiq = fsq = fb = fqg = ^ in Eq- (14), changing the value of f2q- The dependence of the nucleon 
mass on /2q, somewhat "partial derivative" with respect to {0\qq\0) is shown in Fig. 4. The 
physical solution exist only if > 0.99. 

We investigate dependence on the parameters f^g and f^g in the same way. The physical 
solution disappears if f^g exceeds slightly the value f4g = 1 corresponding to the factorization 
hypothesis. There is no minimum corresponding to a physical solution for > 1.01 - Fig. 5. 
On the contrary, the physical solution is not consistent with small values of feg. It requires 



It is instructive also to investigate the dependence of the nucleon mass on the quark scalar 
condensate (0|gg|0) assuming the factorization hypothesis. In this case f^g = /l^, fgg = f2q 
and feq — fig while = 1 in Eq. (14). The result is shown in Fig. 7. There is no minimum 
corresponding to a physical solution since it vanishes for f2g > 1.35. 

3 Inclusion of the radiative corrections of the order cxg 

Here we include the corrections of the order and In in the lowest order of perturbation 
theory. This modifies the contributions to n*(g^). For the contributions of the free quark loop 
Aq and for those of the scalar quark condensate and the four-quark condensate Aq we have 
now [27, 28] 



with Q'^ — —q^. Corrections to the other OPE terms are not included because of the large 
uncertainties of the values of the condensate. 

The numerically large coefficient of the correction of the order ag to the term Aq caused 
doubts in convergence of OPE series [29]. In [28] Eqs. (23) and (24) were used for determination 
of nucleon parameters in framework of the finite energy sum rules. Inclusion of the radiative 
correction was shown to diminish the value of the nucleon mass, assuming that the threshold 
value W"^ does not change. Similar result was obtained in [30] in framework of the Borel 



Ug > 0.96 - Fig. 6. 




(23) 



(24) 
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transformed sum rules. However the authors of [22] demonstrated that the radiative corrections 
modify mostly the values of and W'^, without important influence on the value of nucleon 
mass. 



For the contributions to the Borel transformed sum rules we can write [22] 




( 



M2 ) 




) 



) 



(25) 



with £{x) = 



n-nl 



n=l 



It was shown in [22] that the results do not change much if we put ^^(M^) = ^^(l GeV^). As 
we said above, we assume a;s(lGeV^) = 0.37, which is consistent with the recent data presented 
in Fig. 9. 2 of [26]. A somewhat larger value is given in [21]. Anyway, the nucleon parameteres 
depend weakly on the actual value of ccs-see Fig.2 of [22]. Hence the same referes to the present 
analysis. 

The physical solution is now 



with the contribution of the pole 3 times smaller than that of the continuum. Hence the 

radiative corrections made the unphysical solution even more " unphysical" . Dependence of the 
nucleon mass on the values of condensates is shown in Figs. 2a, 4-7. One can see that the 
physical solution exists in the larger interval of the values of fi than in the case when only 
anomalous dimensions are included. The physical solution (described by Eq. (26) for fb = f^g — 
feq — 1) provides absolute minimum of the function for fb < 1-8, fiq < 1.35 and f^q > 0.2 
- see Figs. 2, 5-7. Dependence on the value f2q without and with factorization assumptions is 
shown in Figs. 4 and 7 correspondingly. 

Note also that the variation of the limits of the Borel window changes the values of nucleon 
parameters only by several percent. For example, solving Eq.(5) in the intervals 0.7 GeV^ < 
M2 < 1, 8GeV 2 we find m = 0.946GeV, while in the interval 1.0 GeV^ < < 1, 6GeV ^ we 
obtain m — 0.953 GeV for = f2q — fiq — 1- There is also an unphysical solution, but the 
corresponding values of are at least ten times larger then for the physical one. 



m = 0.94 GeV, = 2.00 GeV^ W'^ = 1.90 GeV^ . 



(26) 



The unphysical solution is 



m = 0.60 GeV, A^ = 0.56 GeV^, = 0.70 GeV^, 
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4 Summary 



We analyzed dependence of the nucleon mass on the values of QCD condensates in framework 
of the QCD sum rules. We investigated also the dependence of the residue of the nucleon pole 
and of the continuum threshold in "pole + continuum" model of the hadron spectrum. These 
dependences were studied with inclusion of radiative corrections of the order ag and ag \nq^. 

We presented three sets of the results for minimization of defined by Eq. (6) in the Borel 
window defined by Eq.(9). They correspond to total neglection of the radiative corrections, to 
inclusion of the corrections [aglnq^)"^ in all orders and to perturbative inclusion of corrections 
~ as and ccg Incf. 

It is shown that even at relatively small deviations of QCD condensates from the standard 
values the QCD sum rules have an unphysical solution with the contribution of continuum 
exceeding several times that of the nucleon pole. This contradicts the idea of the "pole + 
continuum" model for the hadron spectrum. 

We showed that (neglecting the radiative corrections) at /b < 1 there is an interplay of the 
physical and unphysical solutions. In this case the two solutions can be separated. If both 
unphysical and physical solutions are connected with sharp minima of the dependence of the 
function on nucleon parameters, they can be separated easily. However in the general case, 
for example, for a wide minimum corresponding to the unphysical solution or for a shallow 
minimum of another solution the interpretation of the latter one may be obscure. 

The strongest limits on the values of condensates corresponding to a physical solution emerge 
if one includes the radiative corrections (aslng^)", the limits are weaker if the radiative cor- 
rections are totally neglected and still weaker if the corrections are included perturbatively - 
Figs. 1,3-7. Thus for consistent calculations it is reasonable either to include the radiative 
corrections perturbatively or just to ignore them. 

We demonstrated that for the physical solution the value of the nucleon mass is less sensitive 
to the exact values of the condensates than those of the residue and of the continuum 
threshold W"^. In other words, the uncertainties in the value of condensates influence mostly 
the magnitudes of and W^. 

We found that the nucleon mass depends mostly on the expectation value of the scalar 
quark operator qq. Dependence on the condensates of higher dimensions is much weaker. The 
dependence on the condensate (0|gg|0) becomes weaker if one assumes factorization hypothesis 
for the four-quark and six-quark condensates. 

In the QCD sum rules the nucleon mass obtained a nonzero value due to exchange by 
noninteracting quark-antiquark pairs between the nucleon current and vacuum. As we have 
seen perturbative inclusion of interactions taking place at the distances of the order 1/M ~ 
1 GeV^^, corresponding to inclusion of the radiative corrections ~ Us make the solution more 
stable. The physical solution exists in a broad interval of the values of (OlgglO) if the radiative 
corrections are included perturbatively. On the other hand the unphysical solution becomes 
" more unphysical" , with a stronger dominance of the contribution of continuum over that of 
the pole. This may be important for the QCD sum rules analysis of the many-quark systems 
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carried out nowadays [24, 25]. 
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5 Figure captions 



Fig.l Dependence of the nucleon mass and of per degree of freedom (assuming 1% error 
bar)on value of the gluon condensate /;,. The sohd hnes correspond to the physical solution. 
The dashed lines correspond to the unphysical solution. For convenience the values of are 
reduced by the factor of 5, itlq = 0.93 GeV^ is the value of m corresponding to physical solution 
for fb = 1. For fb > 1.04 only the latter one exists. 

Fig. 2 Behavior of the function m{M'^) for the physical (solid line) and unphysical (dashed line) 
in the case = 0.6. 

Fig.3 Dependence of the nucleon parameters on gluon condensate. Figs. a,b,c show the mass m, 
residue and threshold W'^ correspondingly. Dashed curves show the case with all the radiative 
corrections neglected, dotted curves correspond to inclusion of the anomalous dimensions. Solid 
curves are for perturbative inclusion of the corrections ~ ag- The horizontal dashes denote the 
transitions from the physical solutions to unphysical ones. 

Fig.4 Dependence of the nucleon mass on the parameter f2q. Only the term in Eqs. (11) 
and (25) is modified. The notations are the same as in Fig. 3. 

Fig. 5 Dependence of the nucleon mass on the parameter f4q. Only the term Aq in Eqs. (11) 
and (25) is modified. The notation are the same as in Fig. 3. 

Fig.6 Dependence of the nucleon mass on the parameter /gq. Only the term Bg in Eq. (11) is 
modified. The notations arc the same as in Fig. 3. 

Fig. 7 Dependence of the nucleon mass on the parameter f2q under the factorization hypothesis 
/4g = /Iq) /eg = fiq- The notations are the same as in Fig. 3. 
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